Dynamic Coupling of Convective Flows and Magnetic Field 

during Flux Emergence 

Fang Fang^, Ward Manchester IV^, William P. Abbett^, Bart van der Holst^ 
^ Department of Atmospheric, Oceanic and Space Sciences, University of Michigan, Ann 

Arbor, MI 48109, USA 
^ Space Sciences Laboratory, University of California, Berkeley, CA 94720, USA 



Received 



accepted 



-2- 



ABSTRACT 



We simulate the buoyant rise of a magnetic flux rope from the solar convection 
zone into the corona to better understand the energetic coupling of the solar 
interior to the corona. The magnetohydrodynamic model addresses the physics 
of radiative cooling, coronal heating and ionization, which allow us to produce 
a more realistic model of the solar atmosphere. The simulation illustrates the 
process by which magnetic flux emerges at the photosphere and coalesces to 
form two large concentrations of opposite polarities. We find that the large- 
scale convective motion in the convection zone is critical to form and maintain 
sunspots, while the horizontal converging flows in the near surface layer prevent 
the concentrated polarities from separating. The foot points of the sunspots in 
the convection zone exhibit a coherent rotation motion, resulting in the increasing 
helicity of the coronal field. Here, the local configuration of the convection causes 
the convergence of opposite polarities of magnetic flux with a shearing flow along 
the polarity inversion line. During the rising of the flux rope, the magnetic energy 
is first injected through the photosphere by the emergence, followed by energy 
transport by horizontal flows, after which the energy is subducted back to the 
convection zone by the submerging flows. 

Subject headings: MHD — Sun: interior — Sun: atmosphere 



1. Introduction 



The solar magnetic field emerges at the surface over a wide range of scales, from 
ubiquitous ephemeral regions as small as 10^^ Maxwells, to active reg i ons as large as 10 ^^ 



Maxwells, which emerge at low to midlatitudes (IHagenaar et al. 



2003 



Parnell et al. 



20091 ) 
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Observations from satellites and ground-based telescopes find that the configuration of 
the magnetic field plays an impor tant role in the solar eruptive events, i.e., Coronal Mass 



Ejections (CME) and solar fiares (jCanfield et al 



19991 ). It is thus of great importance to 



study the building up of the magnetic field in the solar atmosphere, as it rises from the 
convection zone. However, the study of the transport of magnetic fiux and energy has been 
hampere d by the inv i sibilit y of subsurface structure s . Wit h the time-distance helioseismic 



analysis. 



KosovichevI ( 119961 ) and 



Zhao &: KosovichevI (120031 ) provide a view on the horizontal 



and vertical fiow velocities on the subsurface layers under sunspots and identify shear fiows 
and rotation of the sunspots underneath the surface, which may build up the energy and 
helicity in the atmosphere. 

Over the past decades, the development of numerical models have also greatly 
improved our under standing of the dynamics and energetics of magnetic fiux emergence. 



Shibata et al 



( 119891 ) describes a two-dimensional magnetohydrodynamic (MHD) simulation 
on the emergence of a horizontal magnetic fiux rop e frorn the photosphere in t o the 



FanI (120081) and 



Jouve fc Brun 



(120091) 



chromosphere using a two-layered atmosphere, 
carry out sets of anelastic MHD simulations on the buoyant rise of the magnetic fiux 
tube from the base of convection zone to near the top, respectively. In particular. 



FanI (120081 ) shows the rotation of the fiux tube driven b y the L orentz force at the two 



ends while the twisted tube is bent. 



Manchester fc Lowl (120001 ) suggested that shearing 



motion driven by the Lorentz force draws the magnetic field parallel to the Polarity 



Inversion 



Fan 



2001 



ine (PIL), which was dernonstrated in simu l ation s of emerging fiux 



Magara fc Longcope 



2003 



Abbett fc Fishei 



2003 



Manchester et al 



ropes 



2004; 



Archontis fc Torok 



2008 



MacTaggart fc Hood 



20091). 



Magara fc Longcopd ( 120031 ) found 



that during emergence, energy fiux through the photosphere is first dominated by the 
vertical fiows while horizontal fiows dominate the late r phase. The energy tra. nsport of shear 



fiows naturally provides an energy source for CMEs (IManchester 



2007 . 



20081). Simulations 
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have also revealed that sh ear flows driven b y the Lorentz force can pr oduce eruptions 



in both magnetic arcades ( Manchester 



2004 



Archontis &: Torok 



2008 



20031) and eniergiii g: flux ropes (IManchester et al 



MacTaggart fc Hood 



20091 ) providing further evidence of a 



mechanism for CMEs, flares and filament eruptions. 

With the availability of more computational resources, numerical models are able to take 



a realistic convection zone overlaid by a self-consistent upper atmosphere ( 


Vo 


gler et al. 


2005 


Stein & Nordlund 


2006; 


Hansteen et al. 


2006; 


Abbett 


2007; 


Gudiksen et al. 


2011 


) . Such 



realistic simulations reveal a w ealth of complex dyn a mic interactions between the magne tic 



field and the convective flows. 



Cheung et al. 



(120071 ) 



Tortosa-Andreu fc Moreno-Insertis 



( I2OO9I ) study the rise of buoyant magnetic flux tubes from the near-surface convection 



zone into the photosphere and chromosphere and flnd the funda mental role of convect ive 



flows in the emergence of th e mag netic flux at the photosphere. 



Martfnez-Sykora et al 



(12008 



Cheung et al 



( 120071 ): 



2OO9I ) simulate an emerging flux tube, which increases the 



size of the photospheric granules and alters the chromospheric structure. More discussion 



on the interaction betw e en th e magnetic fleld and the convective flows can be found in 



Fan 



(I2OO4I ): 



Nordlund et al. 



too% . 



Radiative MHD simulations on a larger scale (iRempel et al 



2009 



C 



report the formation of the sunspots and active regions. 



Rempel et al 



leung et al 



(120091); 



20101) 



Rempel 



( 120111 ) flnd outflows in the penumbral structure dri ven by the Lo r entz 



Cheung et al. 



brces at the surface 



(I2OIOI ) simulates the 



and by the convective flows in the deeper layers, 
formation of a pair of sunspots in an active region with a magnetic semi-torus advected 
through the bottom of the domain and finds the mass removal in the magnetic flux driven 



by the reconnection, and the mi gration of t 



inclusion of turbulent convection (IFang et al 



le flux due to horizontal flows. With the 



20101 ). it is shown that the transport of energy 
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into the corona becomes even more critical as downdrafts return much of the magnetic 
energy back below the photosphere. This simulation treats the emergence of a weak axial 
flux of 3.3 X 10^^ Mx, which is quickly shredded and dispersed to the intergranular lanes by 
the convection motions. 

These radiative MHD simulations illustrate the importance of turbulent convective 



flows in the emergence o: 
expand upon the work of 



the magne t ic str uctures. In light of these previous results, we 



Fang et al. 



(I2OIOI ) by simulating the emergence of a larger twisted 
magnetic flux rope emerging from deeper in the convection zone, and study the energetics 
during its rise through the turbulent plasma to reside in the hot corona. With this work, 
we strive to more deeply understand the energetic coupling between the convection zone 
and the corona provided by the magnetic field. In the following sections, we present the 
results of our simulation. Section |2] describes the numerical methods and the simulation 
steps. In Section [3], we study the results by analyzing the evolution of the magnetic and 
energetic flux under the impact of the convective flows. Finally, we discuss the conclusions 
of our simulation and its implication on further simulations. 



2. Numerical methods 



For our simulation, the MHD equations a re solved wit 
Solar- wind Roe Upwind Scheme (BATS-R-US; iPowell et al. 



1 the B lock- Adaptive Tree 



(119991 )). Additional energy 



source terms are implemented to include the thermodynamic processes in the solar 
atmosphere, i.e. surface cooling at the phot osphere , radi ative cooling of the corona and 



magnetic flux-related heating in the corona ( lAbbett 



20071 ). We describe the radiative losses 



in the low-density corona as optically thin radiation and artificially extend the cooling 



function to treat the so 
by 



Pevtsov et al 



ar surface temperature as well. The coronal heating is approximated 



(120031 ). which gives an empirical relationship between the heating rate 
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and the unsigned magnetic flux at the photosphere. Our model neglects electron thermal 
conduction along magnetic field lines, since the simulations here focus on the transfer of 
energy and magnetic field from the convection zone to the corona. The application of a 



tabular, non-ideal equation of state takes ac count o 



particles in the dense and hot solar interior ( iRogers 



the ionization and the excitation of 



20001). Horizontal boundary conditions 



are periodic, and the lower boundary condition sets the density and temperature of the same 
values found in the initial state of the atmosphere, while keeping the vertical momentum 
refiective ac ross the boundary . The upper boundary is closed. The details of the model are 
described in 



Fang et al 



mm . 



2.1. Background Atmosphere 

To develop the coronal model, we first extract the averaged state of the solar atmosphere 
from the model of 



Fang et al 



(I2OIOI ). at a depth of 2.5 Mm below the photosphere. In 
the convection zone, the plasma is adiabatic as the heating source is absent and radiative 
cooling is negligible in the optically thick medium. The entropy is then nearly invariant 
throughout the deep convection zone. Here the deep convection zone is defined as deep in 
the simulation domain, which covers the upper 10% (20 Mm) depth of the solar convection 
zone. Under the assumption that the plasma is both isentropic and in hydrodynamic 
equilibrium, we solve the following equations for the structure of the atmosphere in the 
deep convection zone with the extracted values providing the upper boundary condition: 



Vp = pg. 



log{^) = C, 



(2) 
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where p, p and g are the plasma pressure, mass density and gravitational acceleration. C 
is the average entropy value in the convection zone at z = —2.5 Mm. 7 is the adiabatic 
index, which can be obtained from the tabular equation of state using the pressure and 
density values, p and p. For the purpose of our simulation here, the atmosphere is extended 
to 21 Mm below the photosphere by solving Equations [T] and [2l The stratification of the 
density and pressure of the extended atm osphere in the convection zone are found to be in 



agreement with that of 



Stein et al. 



(1201 ih 



This solution is then used as an initial condition for a one-dimensional calculation in 
which we apply the surface cooling. The plasma in the simulation domain then goes through 
an hour of thermal relaxation to form a supraadiabatic background atmosphere, which is 
unstable to convective motion. The next step is to perform the full 3-D simulation with the 
size of the domain set to 30 x 30 x 42 Mm'^, extending 21 Mm below the surface and 21 Mm 
above in the vertical directions, while the horizontal area of the domain is comparable to 
the size of a small active region. In the supraadiabatically stratified atmosphere, convective 
motion starts immediately after a small energy perturbation, and relaxes over a period of 
roughly one turnover time of the convection extending to a depth of 21 Mm. After the 
convective motion has relaxed, we impose an initial vertical magnetic field of 1 G to the 
domain, and turn on the heating and radiati ve cooling term in t he coronal region. With the 



20031), the coronal temperature 



magnetic-flux-related coronal heating term ( iPevtsov et al. 
is heated up to 1 x 10^ K within 2 hours. 

The left panel in Figure [T] shows the averaged vertical stratification of the density and 
temperature in the simulation domain. The density drops by about 5 orders of magnitude 
from the bottom of the domain to the photosphere, and the temperature drops by 2 
orders of magnitude. Above the photosphere, the density drops by another 7 orders of 
magnitude, and the temperature increases to 1 MK within 2 Mm. The strong variation in 
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the atmospheric stratification puts large requirements on spatial resolution and restrictions 
on the explicit time step, given the high temperatures and subsequent sound speeds. 
However, the application of the Adaptive Mesh Refinement in BATSRUS greatly facilitates 
the simulation in saving computational resources while keeping the necessary refinements 
to resolve the local turbulent structures. The grid size in our simulations has three levels 
of refinement with cubic cells of size 37.5 km, 75 km, and 150 km. The grid is refined in 
horizontal layers with higher resolution near the surface where temperature and density 
gradients are largest. 

The right panel in Figure [T] shows the variation of the vertical velocity at z = —3 
Mm. Red lines outline the regions with concentrations of Bz field, which coincides well 
with the down-flowing area. The right panel in Figure [2] illustrates the structure of Uz on a 
vertical cut of the convection zone. Small-scale down-flowing plasma merges and forms the 
large-scale persistent downdrafts in the deep convection zone. The large-scale downdrafts 
are very important in the formation of pores during the building up of the active region, 
discussed in Section 13.11 



2.2. Flux Rope 



The stratification of the atmosphere in the simulation domain is maintained self- 
consistently with the implementation of the thermodynamic processes in the model. The 
superadiabatic stratification pr ovides an unstable ba ckground atmosphere for turbulent 



convective motion. As shown by 



Cheung et al. 



(120101 ). convective motion plays an important 



role in the formation of the sunspots and active region. However, the role of the turbulent 
convection in the emergence of the magnetic flux has not been clearly understood yet, 
particularly at larger scales. The aim of our simulation here is thus to study the transfer 
of the energy and magnetic flux of the rope during its rise from the deep convection zone. 
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and its interaction with the surrounding turbulent medium. So after the generation of the 
convection zone with a hot corona, we hnearly superimpose a flux rope upon the ambient 
magnetic field in the deep convection zone, at z = — 10 Mm, indicated by the dot-dashed 
line in the left panel of Figure [H The initial flux rope is centrally buoyan t and twisted 



along the x— axis, as described in iFanI (I2OOII ) and 
equations: 



Manchester et al. 



(I2OOJ) by the following 



B = B 



amb 



(3) 



p = Po(l-r]e--'/''), 



(4) 



1 2 



Po 



(5) 



P = Po(l-^)- (6) 

Here r is the distance from the axis of the rope, and a = 1 Mm, is the radius of Gaussian 
decay, q = —1.5, is the twisting factor. A = 6 Mm, is the length of the buoyant section, po 
is the thermal pressure without the inserted flux rope. Bo = 50 kG, is the strength of the 
magnetic fleld at the axis of the rope, giving a minimum plasma beta value of 50. B^mb 
is the ambient magnetic fleld. The total axial flux on the cross section of the initial flux 
rope is 1.5 x 10^^ Mx. Figure |2] shows the initial flux rope embedded with the convective 
plasma. The left panel illustrates 3-D isosurfaces of the large scale convective flows around 
the flux rope in the deep convection zone, while the right panel shows the structure of Uz in 
the convection zone on the x — z plane with the flux rope. The complexity of the convective 
flows around the flux rope has a strong influence on the emerging process from the deep 
convection zone with large-scale downflows. 
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Results 



The flux rope is initially embedded in the deep convection zone, where it is buffeted by 
convective down- and up- flows. On the way from z = —10 Mm to the photosphere, the 
turbulent flows of various scales signiflcantly reshape and redistribute the magnetic flux. 
On the photosphere, the magnetic flux flrst emerges as many small-area bipoles, which 
then sort themselves to form two large concentrations of opposite polarities via coalescence. 
Figure E] shows the structure of the small active region at 2; = — 1 Mm. It contains two 
pores of opposite magnetic polarities, with smaller-size granules and cooler temperature. 



both remote obsery ations and previous numerical simulations ( iVrabec 



1974 



revealee 


by 


Zwaan 


1985 



Cheung et al. 



20101 ). However, the mechanism that transfers the magnetic flux from the 
tachocline, which is believed to be the origin of the active region magnetic fleld, to the 
solar surface is still to be determined. To date, global scale simulations of the solar interior 
can not yet resolve individual flux tubes. Here, we examine the evolution of the magnetic 
flux rope in the deep convection zone and the photospheric and coronal response to its 



emergence. 



3.1. Formation of Pores 

The horizontal flux rope ai z = —10 Mm, shown by Figure [21 rises in the central section 
due to the depletion of the density and upwelling convective flows. However, the two ends 
of the central section are embedded in large-scale downflows present in the convection zone 
when the flux rope is initiated at t = : 00 : 00. The downflows are illustrated by the blue 
isosurfaces and color contours in Figure |2l Panel (a) of Figure H] shows the structure of 
on the X — z plane at t = 4 : 50 : 00, when the downflows are still present in the convection 
zone at the two ends of the emerged flux rope. The region of the convective downflows 
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in Panel (a) of Figure H] appears in great accordance with the black and green contour 
lines indicating the polarity of the pores. This consistency between the down flowing and 
magnetic concentrated regions suggests a causal relationship between the formation of the 
magnetic pores and the large-scale downfiows. Figure [5] illustrates the temporal evolution 
of the 3-D magnetic field, colored by the Uz value of the local plasma during the rising of 
the central section of the flux rope. The long-lasting, large-scale downfiow, indicated by 
the blue color, drags down the two endpoints of the rising part and fixes them in the deep 
convection zone, forming an fi-shape emerged flux rope within 2.5 hours. The downfiow 
maintains the pores and prevents the two pores of opposite polarities from separation or 
drifting apart, while the central section of the flux rope emerges and expands in the upper 
domain. Thus with the two foot points deeply embedded and fixed in the convection zone, 
the emerged magnetic flux remains highly concentrated in a relatively small area at the 
photosphere. 

In Panel (b) of Figure [3] and Panel (a) of Figure HI sn aall-scale convective granule s 



in the near surface layers appear inside each of the pores. lEmonet fc Cattaned (120011 ) 
also simulates the diminishing horizontal scale of the granules with increasing vertical 
magnetic field. While the magnetic field modifies the convective granules, the downward 
flows produce a bulb of colder plasma. The pressure imbalance with the surrounding 
material thus causes the flux tube to collapse and increases the strength of the magnetic 



field. The flux tube approaches equilibrium again with higher m agnetic 



the surrounding rn aterial. With the convective collapse process (jParker 



pressure balancing 



1978 



Nagata et al. 



Spruit 



1979: 



20081 ) ■ the magnetic field strength at the surface can be increased up to 4 kG, 
much higher than the equipartition field strength with the kinetic energy of the surrounding 
plasma. These near-surface processes, due to the convective motion, play a very important 
role in reorganizing the magnetic fiux after its emergence at the photosphere. The magnetic 
field of the bipolar structures at the photosphere is intensified by the convective collapse in 
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the near surface layers, which is shown by Panel (c) of Figure |H 



3.2. Rotation of the Pores 

The magnetic flux rope travels through the 10 Mm distance and approaches the 
photosphere after 2.5 hours. Figure [6] shows the evolution of the structure of the Bz field 
at the photosphere, with arrows representing the horizontal velocity field. The magnetic 
field concentrates as narrow bands of bipolar fiuxes, as shown in Panel (a) and (b), in 
the regions between the major pores. The horizontal velocity, represented by the white 
arrows, reveals that the bipolar fiuxes are moving in opposite directions, toward the major 
pores of the same polarity. For example, in Panel (a), the rectangle outlines a region with 
fiux emergence, where the negative fiux is moving in the upper right direction into the 
negative pores and the positive fiux bands are moving in the lower left direction into the 
positive pores. The coalescence of the small-scale fiuxes into the major pores facilitates the 
accumulation of the magnetic fiux on the surfac e and, therefore, the formation of the large 



pores shown in Panel (c) and (d). 



Cheung et al. 



(120101 ) simulates the formation of an active 



region and finds that the counterstreaming motion of opposite polarities is driven by the 
Lorentz force. 

The large pores of negative polarity show a coherent pattern of rotation at the 
photosphere after their formation in Panel (a) of Figure IHl The rotation of the pores persists 
during the emergence and the increase of magnetic fiux at the photosphere. However, 
the positive polarity on the left does not present a complete rotation pattern during the 
emerging phase, but the rotation is interrupted by the horizontal motions of the convective 
fiow. The coherent rotation starts to develop on the positive pore after 5.5 hrs, shown by 
Panel (e) and (f) in Figure |6l Figure [7] illustrates the evolution of and the horizontal 
velocity fields ai z = —3 Mm, in the convection zone. Here, we observe a coherent rotation 
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at this depth on the negative polarity as well. The question is then, to what depth does 
the rotation extend. Thus we examine the structure of Uy on the y = plane during the 
rising of the flux rope in the upper panels in Figure [Hi The reversal of the direction of 
Uy in the right side of the domain corresponds to the projection of the rotation of the 
negative polarity on the y = plane. The coherent rotation starts to extend downward at 
t = 4 : 00 : 00 and approaches the depth of -10 Mm in 21 mins. (see Panel (a)). Panel (c) 
shows a very coherent pattern of rotation on the negative polarity at t = 5 : 13 : 00, while 
on the positive pore on the left, the rotation is not obvious. 

Sunspot rotation has l ong been observed and studied in detail an d has been found 



in association with CMEs (IBrown et al. 



2003 



Kazachenko et al. 



20091 ). The ro tation 



mecha n ism for sunspots found at work in our simulation was first described in 



Parker 



(I1979I ). iLongcope fc WelschI (I2OOOI ) simulates the increase of the hehcity of coronal field due 



to the rotatio n of the photospheric footpoints in a dynamic model. Ideal MHD simulations 



by iFanI (|2009| ) illstrate the rotation driven by torsional Alfven waves and the development 
of sig moid-shaped field lin es du ring the fiux e merge nce from the top layer of the convection 



zone. 



Brown et al. 



mom and iMin fc Chael (|2009|) study the behavior of the rotating 



sunspots in solar active regions and found that the rotation speed varies wi th time, rad ius 



and a n gular spacing, with an inc rease in rotation speed in the penumbra. 



Kosovichev 



(|l996|); 



Zhao &: KosovichevI (120031 ) find twists of the magnetic field and vortical fiows around 



the sunspots in the convection zone underneath a rapidly rotating sunspot area. The fiows 
underneath the negative pore in our simulation also exhibit a rotating pattern extending 10 
Mm down into the convection zone. The question is then, what dynamic mechanism during 
fiux emergence causes the observed variation in the rotation speed. Our simulation enables 
a detailed investigation of how sunspot rotation develops in the complex circumstances of 
the convection zone. 
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The development of the coherent rotation is accompanied by the presence of a strong 
Lorentz force in the azimuthal direction (see the upper panels in Figure E]), which is defined 
as fy = jzBx — jxBz- The black and white contour lines in the upper panels of Figure E] 
represent areas with strong negative and positive fy respectively. We find that the Lorentz 
force is driving the coherent rotation as a torsional Alfven wave. At time t = 4 : 21 : 00, 
shown by Panel (a), the Lorentz force is in the same direction as Uy, thus acting to accelerate 
the rotation. However, the Lorentz force reverses direction at time t = 4 : 37 : 00 (see Panel 
(b)) and decelerates the rotation. Panel (c) reveals the structure of the Lorentz force at 
time t = 05 : 13 : 00, which runs in the opposite direction of the rotation. The flux rope has 
rotated past its equilibrium point, and the Lorentz force reverses with the gradient of the 
azimuthal {By) component of the field as the flux rope further emerges. The rotation of 
the pores untwists the field in the convection zone and produces a twisted magnetic field in 
the corona, as shown by the lower panels of Figure El which illustrates the By field on the 
?/ = plane. As such, the rotation also provides a mechanism to build up the helicity and 
magnetic energy in the coronal magnetic field. The white lines in the lower panels outlines 
where the Alfvenic Mach number u/ua = 1- The consistency between the structures of 
the Alfvenic Mach number and the Lorentz force fy suggests that in the pore region the 
magnetic field dominates over the plasma motion, and is responsible for the rotating flows. 

In Figures |6] and [71 the positive polarity does not present a coherent rotation as 
the negative polarity. However, in Panel (c) and (d) of Figure [6l we observe a strong 
flux cancellatio n in the region outl ined by the rectangles from time t = 5 : 00 : 00 to 
t = 5 : 35 : 00. Cheung et al.l (Lolo! ) reports that magnetic reconnection takes place within 
the U-loops formed by the convective downfiows, which can remove the unsigned magnetic 
flux on the photosphere. In our simulation, the total amount of the magnetic flux that 
is cancelled is 10% of the total unsigned flux on the photosphere. The flux cancellation 
strongly interrupts the coherency of the rotation on the positive polarity. The horizontal 
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velocity fields exhibit a converging flow at the two opposite polarities of the cancelled 
flux, both at the photosphere and in the convection zone, yielding a high gradient in the 
magnetic field. The close examination of the flux cancellation event will be the topic of a 
future paper. Another interesting feature we observe in our simulation is that the horizontal 
velocity field in Figure [7] shows a converging horizontal flow field around the magnetic flux 
region. This converging flow constrains the total area of the magnetic concentration, and 
thus prevents the magnetic pores from expansion and separation. 

3.3. Energy Fluxes 

The magnetic flux rope approaches the photosphere within 2.5 hours after its 
initialization, and the photospheric magnetic flux, i.e. the emerged flux, reaches its 
maximum value at time t = 5 : 11 : 00. Afterward, the magnetic pores start to decay slowly 
for the rest of the simulation. However, even at its maximum emergence, the magnetic flux 
rope is far from being fully above the photosphere. This is clearly shown by Panel (b) and 
(c) of Figure HJ where the negative polarity is split into two parts, with one part emerging 
and forming pores while the other is sinking to the deep convection zone. The lower right 
panel of Figure shows the 3-D magnetic field lines of the flux rope when it first approaches 
the photosphere. The red contour color over the field lines represents upflows while blue 
represents downflows. Only about half of the original flux rope is rising with upflows, while 
the other half remains almost stationary or sinks in the convection zone. 

It is of great interest to study the mechanism for transporting the initial magnetic flux 
from the deep convection zone into the photosphere, where the magnetic pores appear, 
and further into the coronal region. Here we calculate the temporal evolution of the total 
unsigned magnetic flux at the photosphere [z = 0), z = —3 Mm and z = 3 Mm, shown by 
the upper panel of Figure [TOl At time t = 4 : 49 : 00, the magnetic flux at z = —3 Mm 
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reaches its maximum value of 2.63 x 10^^ Mx, which is most of the total unsigned flux 
(3.04 X 10^^ Mx) when the flux rope is initially bent. 

At the photosphere, at t = 5 : 11 : 00, the unsigned magnetic flux reaches its maximum 
of 1.37 X 10^^ Mx, 45% of the total initial flux. The unsigned flux at 2; = 3 Mm in the 
lower corona maximizes at t = 05 : 29 : 00 with a value of 8.00 x 10^° Mx. The ratio of 
the emerged unsigned flux with respect to the total initial flux decreases with increasing 
altitudes, from 86% at z = —3 Mm to 26% at z = 3 Mm. The sharp decrease from the 
convection zone to the photosphere is caused by the down flows in the near-surface layers 
that tend to return flux to the deep convection zone. 

The question then, is how the horizontal and vertical flows (driven both by convective 
motions and the Lorentz force) affect the emergence of magnetic flux and the transport of 
the magnetic energy. To address this question we next calculate the magnetic energy flux 
(Poynting flux) passing through three layers of the atmosphere: z = -3, 0, and 3 Mm using 
the following equations: 



shear 



- y ^ (B^u^ + ByUy) B^dS, (7) 

-E^emerge = j ^ ("^^ ~'~ '^z^^' (S) 

The lower panel of Figure Uni illustrates the temporal evolution of the Poynting energy flux 
associated with the horizontal and vertical motions. We find that the energy flux associated 
with the vertical flows at ^ = -3 and Mm, remains negative during the rising and decaying 
phase of the flux emergence. There are several transient positive pulses of energy fluxes 
by the vertical motions dX z = —3 Mm, at times t = 2 : 09 : 00 and t = 2 : 57 : 00. 
These transient positive pulses represent the energy transport when the magnetic flux first 
emerges to the surface with upflows, as shown by the lower right panel of Figure [H Each of 
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these pulses is followed by a sharp increase in the energy flux associated with the horizontal 
flows and a reversal of the energy flux by vertical flows. This time evolution of energy 
shows a process of magnetic flux emergence at the near surface layers with convective 
flows: magnetic flux emerges at the surface as bipoles with upflowing motion, then they 
are quickly pulled apart by the horizontal flows, and concentrate in the downdrafts. This 
process then leads to a positive energy flux by the upflows, followed by a increasing energy 
flux by the horizontal motions, and a negative energy flux associated with the downdrafts. 

Panel (a) and (b) in Figure M show the held at z = -3 and Mm, overlaid by 
white lines showing downdrafts. The concentration of the magnetic flux in the downdrafts 
explains the negative energy flux associated with the vertical flows in the lower panel from 
3.0 hrs to 5.0 hrs while the total magnetic flux is still increasing in the upper panel. During 
this time period, more than 60% of the total unsigned flux is in the downflowing region for 
both z = - 3 and Mm layers. At time t = 5 : 11 : 00, when the photospheric unsigned flux 
maximizes, 70% of emerged flux is concentrated in the downdrafts. Panel (c) and (d) show 
the Poynting fluxes at the photosphere associated with the vertical and horizontal flows, 
deflned as: 

-^shear = — {BxUx + ByUy) B (9) 

47r 

^emerge = ^ [BI + B^) U,. (10) 

Femerge IS negative in the magnetic polarities while postive in the areas between the 
polarities, where the magnetic flux is emerging. In the pores, the energy flow is instead 
dominated by the Fghear, shown by Panel (d). At 2 = 3 Mm, in the corona, the energy 
flux by the vertical flows shows short periods of postive values. The increased energy 
flux is driven upward by convectively-driven magnetoacoustic sh ocks, and could be 



interpreted in terms o 



Martfnez-Sykora et al 



;he d ynamics of Type I spicules (see e.g., iHansteen et al.l (120061 ): 



f l2009h ) 
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The Poynting flux associated with the horizontal flows dominates the energy transport 
during the flux emergence. At the photosphere, the total energy transport by this flux is 
1.35 X 10^^ ergs, while vertical flows transport 5.77 x 10'^^ ergs, 42% of the energy back to 
the convection zone by the end of the simulation. The horizontal flows here include the 
rotation of the magnetic pores, the separating motion of the small bipoles and the shearing 
motion along PILs. The rotation of the pores is discussed in Section 13. 2[ which transports 
both the magnetic energy and helicity into the corona region. The extension of the rotation 
motion into the deep convection zone shown by Figure [8] greatly impacts the spatial 
distribution of the Poynting flux by the rotation. Figure E] shows the separating process of 
the small bipoles, which tend to move apart from each other after their emergence at the 
photosphere and merge into large pores with the same polarity via this self-sorting process. 
The horizontal separating flow on the small bipoles builds up energy in the near surface 
region. In Panel (c) and (d) of Figures |6] and [TJ the black rectangles outline the area with 
a large-scale magnetic flux cancellation. The converging flows across the PIL build up the 
magnetic gradient in this area, and the shearing flow along the PIL, shown in Panel (c) and 
(d) of Figure El drives the magnetic field lines along the PIL and builds up a highly sheared 
magnetic field configuration, w hich plays an important role in eruptive e vents such as fla res, 



filament eruptions and CMEs ( iMikic &: Linker 



1994 



Wu fc Guo 



1997 



Antiochos et al. 



19991 ). A detailed analysis of this flux cancellation event will be the topic of a future paper. 



4. Discussion and Conclusions 

In Section 13.11 and 13.21 we examined the structures of the vertical and horizontal 
velocity flelds separately and their influence on the emerging magnetic flux rope. The 
persistent large-scale downflows help form and maintain the bipolar pores of uniform 
polarity in the deep convection zone, while the small-scale downflows in the near-surface 
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lay ers inte n sify t he strength of the emerged flux tube by convective collaps e , as s uggested 



by 



Parkerl (119781 ) and found in previous simulations by 



Stein fc Nordlundl fl2006h . At 



the photosphere, the horizontal flows act on the small newly-emerged bipoles, to not 
only separate opposite polarities, but also sort them into large pores of uniform polarity. 
This coalescence of small-scale fluxes enables the formation of the magnetic pores with 
a unsigned flux of up to 1.37 x 10^^ Mx, which is limited by the scale of our simulation. 
The horizontal flows of the magnetic pores exhibit a coherent rotation extending down 



to the deep convect ion zone, which is c 



earlier simulations f lLongcope fc Welsch 



early driven by the Lorentz force as was found in 



2000 



Fan 



20091 ). Our simulation illustrates that 



this rotation mechanism operates in a realistic convection zone but requires the coalescence 
of a well formed pore and that complex interactions that break the symmetry of the pore, 
(such as flux cancellation) may disrupt the rotation. 



The energy transport due to the horizontal an d vertical flows is calcu. 



Abbett fc Fisheil (120111 ) finds the 



ated in the 



domain from the convection zone into the corona. 
Poynting energy flux flows into the interior below the visible surface and flows into the 
corona above it and suggests surface convection as the energy source of the separatrix in the 
energy flux. Our study shows a negative total energy flux aX z = —3 Mm, and positive at 
z = and 3 Mm, which agrees with this result. The energy flux is initially dominated by the 
emerging flow on the magnetic flux, quickly after which opposite polarities separate, giving 
rise to an increase in the energy flu x associated with the horizo n tal flows. This geri e ral tr end 



has been found in early work e.g. 



Fang et al 



Magara fc Longcopd (120031 ): 



Manchester et al. 



(I2004h : 



(120101 ). but these works did not illustrate the extent to which downflows control 
energy transport in the convection zone. In the convection zone and at the photosphere, 
the flux concentrates in the downflow drafts, yielding a negative energy flux associated 
with the vertical motion. The horizontal flows are thus the dominating carriers of the 
energy transport into the corona. We identify three types of horizontal flows in Section 13.21 
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the rotation, the separating motion of the bipoles, and the shearing motion along PILs. 
The rotation and shearing both transfer the magnetic energy and hehcity into the corona. 
And it is the separating and shearing motions that build up the magnetic energy in the 
near-surface region, which is crucial for the eruptive events. 

In our simulation, we also find an area with cancellation of magnetic flux up to 10^° 
Mx within 0.5 hour, shown by Panels (c) and ( d) of Figure El The coronal response to this 



Green et al 



(I2OIII ) finds the evolution of 



large-scale flux cancellation is of great interest, 
coronal fields into a highly sheared arcade and then a sigmoid in an active region where 1/3 
of the flux is cancelled. Whether the cancelled flux in our simulation produces a flux rope 
in the coronal region will be the topic of a future paper. Another interesting feature is the 
dec ay of the magnetic pores from tim e t = 6 : 00 : 00 to t = 8 : 00 : 00, shown by Figure 



[TUl Ivan Ballegooijen fc Mackayl (j2007l ) suggests that the submerged magnetic field repairs 
the toroidal flux ropes from which the initial flux emerged. The decayed magnetic flux thus 
may play a fundamental role in replacing the magnetic field in the solar interior. 
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Fig. 1. — Left: The vertical stratification of the density (sohd) and the temperature (dashed) 
of the solar atmosphere. The dash-dotted line indicates the initial location of the inserted 
flux rope. Right: The vertical velocity structure ai z = —3 Mm. Red lines show regions 
with \Bz\ greater than 10 G. 




Fig. 2. — The initial structure of magnetic flux rope with color showing the vertical velocity 
in the convection zone. Left: magnetic flux rope embedded in the convective flows with 
blue isosurfaces indicating upflows at 1 km s^^ and red downflows at -1 km s~^. Right: 
X — z plane of the convection zone with the inserted flux rope. 
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Fig. 3. — The structure of Bz (a), Uz (b) and T (c) at z = —1 Mm plane at t = 5:09:00. 
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Fig. 4. — The structure of Uz (a), Bz (b) and Bz (c) in the y = plane at t = 4:50:00 (a 
and b) and t = 5:22:00 (c). The black and green lines in Panel (a) outline Bz = -1, 1 kG, 
respectively. The dash-dotted line in Panel (b) and (c) indicate the initial location of the 
axis of the inserted flux rope. And the black lines in Panel (b) and (c) indicate the magnetic 
field lines by ignoring the By component to show the direction of the field confined to the 
X — z plane. 



Fig. 5. — The temporal evolution of 3-D magnetic field lines colored by local Uz values with 
red indicating downflows and blue upflows from -2 km s~^ to 1 km s~^, at t = 1:05:00, 
1:29:00, 1:49:00 and 2:29:00. 
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Fig. 6.— The structure of B, field at z = Mm at t = 3:45:00 (a), 4:15:00 (b), 5:10:00 (c), 
5:35:00 (d), 6:23:00 (e) and 7:41:00 (f). The arrows show the horizontal velocity field. 
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Fig. 7.— The structure of B, field at z = -3 Mm at t = 4:40:00 (a), 5:00:00 (b), 5:10:00 
(c), and 5:30:00 (d). The arrows show the horizontal velocity field. 
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Fig. 8. — The structure of Uy (upper panels) and By (lower panels) at time t = 4:21:00 (a 
and d), 4:37:00 (b end e) and 5:13:00 (c and f) in the y = plane. The black and white 
contour lines in the upper panels represent regions with Lorentz force density fy = -20 and 
20 dyne cni'^ respectively. And the black lines with arrows in the lower panels indicate 
streamlines of [B^^Bz) while white lines outline Alfvenic Mach number u/ua = 1- 
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Fig. 9.— The structure of B, field at z = -3 Mm at t = 3:40:00 (a), B, field at z = 
Mm at t = 4:03:00 (b) and Poynting fiuxes associated with vertical (c) and horizontal fiows 
(d) aX z = Mm at t = 4:03:00. The white lines in Panel (a) and (b) indicate downfiowing 
regions with Uz = -0.1 km s~^, and the black lines in Panel (c) and (d) show regions with 
\Bz\ = 500 G. 
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Fig. 10. — The temporal evolution of magnetic flux (upper) and energy flux (lower) asso- 
ciated with the horizontal flows (solid) and vertical flow (dashed) at z = —3 (purple), 
(green) and 3 (red) Mm from time t = 2 : 00 : 00 to t = 8 : 00 : 00. 



